Phases of Dense Quarks at Large N c 
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Abstract 

In the limit of a large number of colors, N c , we suggest that gauge theories can 
exhibit several distinct phases at nonzero temperature and quark density. Two are 
familiar: a cold, dilute phase of confined hadrons, where the pressure is ~ 1, and a 
hot phase of deconfined quarks and gluons, with pressure ~ N~. When the quark 

■ chemical potential /j, ~ 1, the deconfining transition temperature, Td, is independent 
. of fj,. For T < Trf, as fi increases above the mass threshold, baryons quickly form 

a dense phase where the pressure is ~ N c . As illustrated by a Skyrme crystal, 
chiral symmetry can be both spontaneously broken, and then restored, in the dense 
phase. While the pressure is ~ N c , like that of (non-ideal) quarks, the dense phase 
Q\ | is still confined, with interactions near the Fermi surface those of baryons, and not 

of quarks. Thus in the chirally symmetric region, baryons near the Fermi surface 
are parity doubled. We suggest possible implications for the phase diagram of QCD. 
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O ■ 1 Introduction 



Many of the observed properties of QCD can be understood, at least qual- 



itatively, by generalizing from three to a large number of colors, N c — > oo. 
[1,2,3]. For example, consider the phase transition at a nonzero temperature, 
T [4]. At low temperature confinement implies that all states are color singlets, 
such as mesons and glueballs, with a pressure ~ N® ~ 1. At high tempera- 
ture gluons in the adjoint representation deconfine, contributing ~ N% to the 
pressure. Thus one can define the deconfining transition simply by the point 
where the term ~ in the pressure turns on. The transition temperature 
for deconfinement, T^, is expected to be of order one at large N c , on the order 
of a typical QCD mass scale, such as the renormalization mass parameter, 
A.qcd ~ 200 MeV. Arguments suggest that deconfinement is a strongly first 
order transition, with a latent heat ~ N%. Since the free energy of Nf flavors 
of deconfined quarks is ~ N c Nf in the limit of large N c , and small Nf, decon- 
finement probably drives chiral symmetry restoration at T^. Several of these 
features have been confirmed by numerical simulations on the lattice [5]. 
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In this paper we consider the phase diagram in the plane of temperature 
and quark chemical potential, \i [6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24]. 
It is usually assumed that for all T and /x, there is a single transition, at which 
both deconfmement and chiral symmetry restoration occurs, shaped something 
like a semi-circle. If the transition is crossover for \x = and T ^ [17], there 
could be a chiral critical end point in the plane of T and /i [13]. From numer- 
ical simulations on the lattice [17,18] within errors the two transitions always 
appear to coincide, at least for T^O and small fi. 

At large N c , we find a very different phase diagram in the T-/i plane, 
in which the deconfining and chiral transitions split from one another. For 
/j, ~ 1, the deconfining transition temperature is independent of \i. At low 
temperatures and densities, there is the usual confined phase of hadrons, with 
chiral symmetry breaking. The confined phase is baryon free, as a Fermi sea of 
baryons first forms at a value near the lightest baryon mass. Within a narrow 
window in fi, ~ l/-^c> there is then a rapid transition to a dense phase, with 
a pressure ~ N c . The properties of dense phase(s) are illustrated by a Skyrme 
crystal [25,26]. Although the total pressure is ~ N c , like that of quarks, the 
dense phases are confined, with interactions near the Fermi surface dominated 
by baryons. We suggest that the dense phase undergoes a chiral transition for 
/i ~ 1. In the chirally symmetric phase, the baryons are parity doubled [27], 
consistent with the constraint of anomalies at nonzero density [28]. 

Admittedly, all of our arguments are merely qualitative. Even so, we think 
it worth pursuing them, because they are so different from naive expectation. 
Of course our analysis could simply be an artifact of the large N c expansion, 
and of limited relevance to QCD, where N c = 3. At the end, we suggest what 
our analysis might imply about the phase diagram of QCD. 

2 Review of Large N c 

If g is the gauge coupling, the 't Hooft limit is to take N c — > oo, holding g 2 N c 
fixed [1,2]. This selects all planar diagrams of gluons. Holding Nf fixed as 
N c — > oo, quark loops are suppressed, and the only states which survive have 
a definite number of quarks and anti-quarks. 

Mesons are composed of one quark anti-quark pair, and are free at infinite 
N c : cubic interactions vanish ~ l/^/Wl, quartic interactions, ~ l/iV c , etc. 
Glueballs are pure glue states, with no quarks or anti-quarks; their cubic 
interaction vanish ~ 1/jV c , and so on. Except for Goldstone bosons, the 
lightest bosons have masses ~ h-QCD- 

In contrast, baryons are rather nontrivial. To form a color singlet, they 
have N c quarks. Assuming that each quark has an energy of order Aqcd, 
the mass of a baryon M B ~ N c Aq C d [2,3]. For instance, a gluon exchanged 
between any two quarks contributes ~ g 2 N 2 ~ N c , and contributes to an 
average Hartree potential. Care must be taken with diagrams to higher order: 
one includes diagrams which modify the average potential, but not iterations 
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thereof. Thus two gluons exchanged between two different pairs of quarks is 
~ N 2 , but represents an iteration of the average potential. A diagram which 
modifies the average potential is given by three quarks which emit three gluons, 
and which then interact through a three gluon coupling, fig. (41) of [2]; this 
is ~ g A Nl ~ N c . 

The scattering of two baryons begins with the exchange of two quarks 
between the baryons, with a gluon emitted between them, fig. (35) of [2]. 
Since each quark can have a different color, this diagram is ~ g 2 N 2 ~ N c . It 
is also possible to exchange quarks of the same color without gluon emission; 
this is also ~ N c , fig. (36) of [2]. As for the average potential in one baryon, 
the two body scattering amplitude remains ~ N c to higher order in g 2 . To see 
this, it is necessary to pick out contributions which are two baryon irreducible, 
from those which represent iterations of the two baryon potential. 

Continuing, the scattering between three baryons is ~ N c : if a quark 
in each baryon emits a single gluon, which interact through a three gluon 
interaction, the amplitude is ~ g A N^ ~ N c . Again, multiple gluon exchange 
does not produce higher powers of N c , once one picks out interactions which 
are three baryon irreducible. 

In general, the scattering of Ai baryons is of order ~ N c . Thus unlike 
mesons or glueballs, whose interactions vanish at infinite N c , baryons interact 
strongly, with couplings of strength ~ N c . 

Of course baryon interactions can also be viewed as arising from the ex- 
change of color singlet mesons. The coupling between a meson and a baryon 
~ y/N~ c , and so single meson exchange gives a two baryon interaction ~ N c , 
as above. While it appears that multiple meson exchange will lead to higher 
powers of N c , this does not occur, due to an an extended SU(2Nf) symme- 
try [3]. These cancellations are subtle, and surely have analogies in nuclear 
matter. For our purposes, however, all we require is that the scattering of M. 
baryons is always of order ~ N c . 

The SU(2Nf) symmetry implies the low energy spectrum of baryons is 
highly degenerate. The lowest mass baryons form multiplets of isospin, /, 
and spin, J [3]. These multiplets have I = J, from 1/2 to N c /2 for odd N c . 
(For QCD, there is one such state, the A.) The splitting in energy between 
the states in these multiplets is of order Mb ~ M(l + kJ 2 /N 2 ), where k is a 
constant. These are the lightest states: there are other excited baryons with 
masses ~ Aq C d above the lightest. 

At zero temperature, there is no Fermi sea until the chemical potential 
exceeds the mass of the fermion. Let M be the mass of the lightest baryon, 
M ~ N c [2]. It is natural to define a "constituent" quark mass, 

m q = M/N c , 

which is of order one at large N c . Thus at T = 0, there is no Fermi sea until 
the baryon chemical potential fis > M; for the quark chemical potential, this 
is ji > m q . 
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To illustrate how quarks enter at large N c , consider the gluon self energy 
at nonzero T and fi. To lowest order in g 2 , at zero momentum this is gauge 
independent, equal to the square of the Debye mass. For Nf massless flavors, 
its trace equals 



IF"(0) = g z 



c+ 2 ) 3 + 2vr2 



(1) 



Taking N c — - > oo, holding g 2 N c fixed, we see that the gluon contribution, 
~ g 2 N c T 2 ~ T 2 , survives. This is the first in an infinite series of planar, gluon 
diagrams at infinite N c . In contrast, whether for T^O and fi ^ 0, the quark 
contribution is only ~ g 2 , and so suppressed by ~ 1/N C . 

This is true order by order in perturbation theory, both in vacuum and 
for all T and fi ~ 1: holding Nf fixed as N c — > oo, the effects of quarks 
loops are suppressed by ~ 1/N C [1,2]. This is simply because there are ~ N 2 
gluons in the adjoint representation, but only ~ N c quarks in the fundamental 
representation. Since the quark contribution, relative to that of gluons, is 
~ Nf/N c , it is essential to hold Nf fixed as N c — > oo; i.e., to take of limit of 
large N c , but small Nf. 

In this limit, we can immediately make some broad conclusions about the 
phase diagram in the T — fi plane. At fi = 0, one expects that the deconfining 
transition temperature T d ~ Aq C d [4], which appears to be confirmed by 
numerical simulations on the lattice [5]. Since quarks don't affect the gluons, 
the deconfining transition temperature is then independent of fi, Td(fi) = Td(0) 
for values of fi ~ 1. This is illustrated in fig. (1): in the plane of T and fi, 
the phase boundary for deconfinement is a straight line. The theory is in a 
deconfined phase when T > Td, and in a confined phase for T < T^. 
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Fig. 1. Phase diagram at infinite iV c in the plane of temperature and quark chemical 
potential. The blue line in the quarkyonic phase indicates a guess for the position 
of the chiral phase transition. 

In fact, consider the "box" in the lower, left hand corner of the T — fi 
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plane, where T < Tj and /x < m q , fig. (1). For particles of finite mass, as long 
as T ^ 0, one expects some population of fermions in the thermal ensemble. 
At large N c , however, the baryons all have a mass ~ N c ; thus if /x < m q , as 
long as T < the relative abundance of baryons is ~ exp(— kN c ), with k a 
number of order one, and so the baryon abundance is exponentially small at 
large iV c . (There is a small window in which baryons can be excited, when 
xa — m q ~ l/iV c , where it costs ~ 1, and not ~ N c , to excite baryons.) 

This box, T < T d and /x < m q , is the usual, confined phase of hadrons. 
At nonzero temperature, the pressure ~ 1 is due exclusively to mesons and 
glueballs, with only exponentially small contributions from baryons. Even for 
/x > m q , the only baryons are those in the Fermi sea, or excitations thereof. 
Simply because they are too heavy, virtual baryon anti-baryon pairs never 
contribute at large N c . This is not true if /x grows like a power of N c , but we 
generally do not consider this regime, except following eq. (9). 

Henceforth we concentrate on cool, dense quarks: remaining in the confined 
phase, T < T^, and moving out in /x from m q . This is the "quarkyonic" phase 
in fig. (1). We explain this terminology later, but stress that as it occurs 
for T < Td (and \i ~ 1), that it is confined. We note that Cleymans and 
Redlich showed that in QCD, phenomenologically the boundary for chemical 
equilibriation begins at /x ~ m q = M^/3 = 313 MeV when T = [14], 
reminiscent of fig. (1). 



3 Narrow Window of Dilute Baryons 

We start by working at zero temperature, very close to the point where a Fermi 
sea of baryons forms. The Fermi momentum for baryons, kp, is k 2 F + M 2 = fi 2 B . 
If kp is (arbitrarily) small, we have an ideal gas of baryons. For such a gas, 
the baryon density is n(kp) ~ k F , the energy is e ~ k F /2M, and the pressure 
is 

k 2 1 k 5 

Pideal baryons ~ 7l{k F ) — ~ — - . (2) 

M l\ c Aqcd 

For such a dilute gas of baryons, the pressure at T = is very small, ~ l/N c . 

Now consider increasing fcp further, so that the additional resonances of 
the baryon condense. There are several effects which enter. The first is to 
include the resonances with I = J, representing the generalization of the A, 
etc., at large N c . Each species contributes to the pressure as (kp — (k// 2 + 
kjJ 2 )Aq CD ) 5 / 2 /M. Summing over spin and isospin gives a total contribution 
of order 

Thus while the sum over resonances changes the dependence upon kp, it still 
contributes to the pressure ~ 1/N C . 

In contrast, once the nucleon Fermi momentum increases, the effect of 
interactions quickly dominates. The amplitude for the four point interaction 



5 



between baryons includes a term ~ Nc^ip) 2 / Aq Cd . The coupling for this 
interaction has dimensions of inverse mass squared, which we assume is typical 
of Aqcd- At low densities this interaction contributes of order density squared, 
or 

Ap „ <k F f k% 

O-ftwo body int.'s ~ ^ c ~ ^ c ■ \V 

QCD A QCD 

Likewise, six point nucleon interactions contribute as density cubed, or ~ 
N c k 9 F /A.Q CD , etc. 

Now clearly one cannot trust this series when the nucleon Fermi momentum 
is of order Aqcd- What is interesting is when the series breaks down. Consider 
balancing (2) and (4): the two terms are comparable when kp/N c ~ N c k F , or 

c 

Thus at very low nucleon Fermi momentum, kp ~ 1/-^, two body nucleon 
interactions are as important as the kinetic terms. Contributions from reso- 
nances, (3), are suppressed by one factor of the density, kp ~ 
three body interactions. 

As kp increases beyond this point, contributions from the Fermi sea, even 
including the increasing number of states, are irrelevant. Instead, when kp ~ 
Aqcd, the pressure is completely dominated by baryon interactions. Since the 
interactions between M. baryons are ~ N c , all baryon interactions contribute 
equally, to give a dense phase in which the pressure is ~ N c . 

It is important to stress that at large N c , the window in which baryons are 
dilute is very narrow. We expect to enter a dense phase, with pressure ~ N c , 
when the baryon Fermi momentum kp ~ 1. In terms of the quark chemical 
potential, though, 

fJ-B-M fc| 1 2 
u — m a = = ~ — - k p . o 

P q N c 2MN C iV c 2 F K ' 

That is, for /i one enters a dense regime within 1/Nl of the mass threshold. 
This is why in fig. (1), we have indicated that the quarkyonic phase begins 
right at m q : the window in which one has dilute baryons is only ~ l/A^ in 
width. 

This discussion is somewhat naive. If the potential between two nucleons 
is attractive at large N c — as it is in QCD — then at arbitrarily low densities 
free nucleons collapse to form bound nuclear matter. (This is analogous of the 
tendency of nuclear matter in QCD to go to the most stable state, which is 
one of iron nuclei.) The mass threshold for baryons is then not at = Mb, 
but at — Mb — 8E, where 5E is the binding energy of nuclear matter at 
large N c [29]. One expects 5E = N c 5e, with 5e ~ Aq C d- Then in fig. (1), the 
quarkyonic phase begins not at m q , but within ~ l/N% of \i = m q — be. 

In QCD, the nuclear binding energy is anomalously small, 5e = 5E/3 ~ 
5 MeV, versus m q = 313 MeV. We do not know if 5e is generically small in 
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the limit of large N c . If so, it is surely related to the SU(2Nj) symmetry [3]. 
In QCD, however, this may be an "accident" of 2 + 1 light flavors. 

In QCD, a gas of nucleons doesn't directly collapse to nuclear matter, 
but instead exhibits a liquid-gas transition. This is because of two effects, 
the kinetic energy of individual nucleons and Coulomb repulsion. Neither is 
important at large N c . As seen above, effects of kinetic energy are automat- 
ically ~ 1/N C . Likewise, taking the baryon electric charge to be of order one 
at large N c (as is necessary for electromagnetism to remain weakly coupled), 
then Coulomb repulsion is ~ 1, and negligible relative to the nuclear potential, 
~iV c . 

4 Dense Baryons as a Skyrme Crystal 

In the previous section we saw that when the nucleon Fermi momentum 
k F ~ Aqcd, that one goes into a dense phase, dominated by baryon-baryon 
interactions. To understand what might happen in this regime, in this section 
we review the properties of Skyrme crystals [25,26]. Our principal interest is 
as an example of a confined theory which nevertheless has a chirally symmetric 
phase. 

The usual Skyrme model is a sum of two terms, 

L = /^tr|V^| 2 + Ktr[V^K] 2 , V„ = tfd„U , U = exp(i 7 r// 7r ) , (7) 

where f n is the pion decay constant, and k a coupling constant, with f% ~ n ~ 
N c . We limit ourselves here to the Skyrme model for two flavors, where tt is the 
pion field. There are many other terms besides those in (7); terms from the 
anomaly contribute through the Wess-Zumino-Witten (WZW) Lagrangian. 
The two terms above must be viewed as the leading terms in a derivative 
expansion. Terms with higher numbers of derivatives have coupling constants 
with dimensions of inverse mass squared. The mass dimension of these other 
terms is presumably set by (inverse) powers of Aq C d- 

A single Skyrmion is given by a solution to the field equations from (7) 
over all of space-time. Since the terms in the action are ~ N c , the energy of 
a configuration is also of the same order, and represents a single baryon. 

As shown by Klebanov [26], a realistic crystal is given by considering peri- 
odic solutions in a finite box. Like the energy of a single baryon, the energy of 
the Skymrion crystal is automatically ~ N c , with one baryon per box. Solving 
the Skyrme equations of motion for a system with cubic symmetry is tech- 
nically involved. For many crystals, however, it is known that a reasonable 
approximation is to chop off the corners of the cube, and to consider the the- 
ory on a sphere. This approximation was adopted by Kutschera, Pethick, and 
Ravenhall, and also by Manton [26] . Many properties of the crystal with cubic 
symmetry are especially transparent for a spherical geometry. 

If R is the size of the sphere, the solution is constructed so that there is 
one baryon per sphere, so the baryon density is l/(47r_R 3 /3). A crystal is a 
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bad approximation for large R, but presumably reasonable when R ~ y/n/ f n . 
At large N c , this mass scale ~ 1. 

For large spheres, the chiral symmetry is broken, as the U field points in 
a given direction in isospin space; typically, n — > at spatial infinity. As R 
decreases, the stationary point is distorted by the finite volume of the sphere. 

As the radius of the sphere becomes small, there is a phase transition to a 
chirally symmetric phase. For small spheres, the stationary point is just the 
identity map, from from S 3 of the U's to S 3 of space, taking ir a ~ r a . This 
has unit baryon number per spherical volume, but it is also easy to see that 
the integral of U, over the sphere, vanishes. 

The restoration of chiral symmetry is less obvious for a crystal with cu- 
bic symmetry. It follows from the half-Skyrmion symmetry of Goldhaber 
and Manton, where the total chiral order parameter cancels between different 
regions of the crystal [26]. (For more than two flavors, the correct represen- 
tation of the chirally symmetric phase is presumably a generalization of this 
half-Skyrmion symmetry.) 

The crucial test for the restoration of chiral symmetry is that the excitation 
modes fall into chiral multiplets of the unbroken symmetry. This was shown 
for the mesonic excitations of the crystal by Forkel et al. [26]. Since the 
baryon current is topological, the baryon excitations are not evident. In a 
chirally symmetric phase, they must be parity doubled. In detail, this happens 
because the baryon is a topological current. Thus it is given by integrating over 
the entire box; if the configuration is chirally symmetric, so are the integrals 
thereof. 

The Skyrme crystal does not give one insight into all properties of the 
system. The pressure of the system is not obvious, nor even the chemical 
potential of baryons. What one can do is to compute how the energy, per 
cell, depends upon the density. If one takes only the two terms in the Skyrme 
Lagrangian of (7), then the term with four derivatives dominates at small 
R, which is high density. Since this term is scale invariant, one automatically 
finds that the relationship between the energy density, e, and the density, n, is 
that for a conformally symmetric theory, e ~ n 4//3 , controlled by the coupling 

H . 

This is an accident of keeping only two terms in the Skyrme lagrangian. 
Terms with six derivatives, for example, are also proportional to N c , with 
dimensions of 1/Aq Cd . When the size of the crystal is ~ 1/Aqcd, however, 
all such interactions are equally important. This is analogous to the count- 
ing for baryon baryon interactions, which are always ~ N c , and which are 
characterized by mass scales ~ Aqcd- 

We stress that the Skyrme crystal is only meant as an illustration. For 
example, while naively one expects that a system of heavy particles forms 
a crystal, since the interactions are as large, also ~ N c , even this may not 
necessarily follow. 
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5 Quarkyonic Matter at Large N c 



Skyrme crystals illustrate what might happen as \x increases from m q (or \xb 
from M). Alternately, one can consider working down in fi. For very large /i, 
fi ^> Aqcd, one should be able to compute the total pressure by perturbation 
theory, in the QCD coupling g 2 : 

P pert .M ~ N c N f /i 4 F (g 2 (»/A QCD ), N f ) , (8) 

The function Fq has been computed to ~ g A [6]. At large N c , it is given by 
planar diagrams with a single quark line, and an arbitrary number of gluon 
insertions. The pressure of ideal (massless) quarks is ~ /i 4 . Perturbative 
corrections are given by a gluon knocking a quark in the Fermi sea out of it, the 
subsequent rescattering of the quark(s) and gluons, etc. When ji ^> Aq C d, for 
quarks deep in the Fermi sea, all such contributions are due to hard scattering, 
with the energy and momenta transferred ~ /i. Since the gluon has momentum 
components of order /j, the the coupling runs according to this scale, with the 
/3-function of the pure glue theory, (8). Thus in perturbation theory, the 
pressure is a power series in ~ 1/ 'log(/j/ Aqcd) > etc. times /i 4 . 

To be able to compute the total pressure reliably in perturbation theory, 
we only need to assume that /i ^> Aqcd- In particular, it is not necessary to 
assume that /i grows like a power of N c . Presumably a perturbative calculation 
is applicable for, e.g., fi > 10 2 Aqcd- When N c is absurdly large, such as 
N c = 10 12 , this is still smaller than any other scale which enters at large N c , 
such as Nc^Aqcd ~ W 3 Aqcd] see the discussion following eq. (9). 

Now consider pushing the perturbative computation of the pressure down 
to /i ~ Aqcd- Non-perturbative contributions to the pressure enter, through 
terms such as ~ /i 2 , eq. (10). Even so, outside of the window of dilute baryons, 
sec. 3, baryons are dense, and we expect that the pressure remains ~ N c . 

This then raises the central conundrum of our work: for T < T^, and \i ~ 1, 
it appears that we can describe the system either as one of confined baryons, 
or as one of quarks. Admittedly, the baryons interact strongly, ~ iV c ; likewise, 
and especially for fi ~ Aqcd-, the Fermi sea of quarks is far from ideal. But 
how can both pictures apply? 

We suggest the following resolution. At large /i 3> Aq C d, for quarks far 
from the Fermi surface, their scattering can be reliably computed in perturba- 
tion theory. This is reasonable: at large /i, the density of quarks per hadronic 
volume ~ 1/Aq CD , is large. In such a dense medium, a quark doesn't know 
which baryon it belongs to; then, for the most part, it is appropriate to view 
the system as one of (non-ideal) quarks. 

Even at large /i, however, it is essential to consider separately the scattering 
of particles within ~ Aqcd of the Fermi surface. In this regime, quarks interact 
by exchanging gluons with momenta ~ Aqcd- At infinite iV c , where quarks 
cannot screen gluons, we know how quarks at fi ^ scatter: exactly as for 
/i = 0. When T < T d , then, the theory is in a confined phase, and so near the 
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Fermi surface, it is appropriate to speak not of the scattering of quarks, but 
of baryons. 

We term this a "quarkyonic" phase: a quark Fermi sea, with a bary- 
onic Fermi surface. The width of the baryonic surface is ~ Aqcd, so when 
yU ~ Aqcd, it is all baryons. As /i increases, the baryons form a band, of 
approximately constant width, on the edge of the Fermi surface. There is 
no quantitative difference between a quarkyonic phase, with a wide baryon 
surface, and one with a narrow surface: they smoothly interpolate from one 
to the other. At large N c , it is possible to differentiate the quarkyonic phase, 
with a pressure ~ N c , from that in the hadronic phase, ~ 1, or the deconfmed 
phase, ~ N%- This clear distinction is only possible at large N c (and small 



The effects of a baryonic Fermi surface show up in the pressure through 
terms which are powers of ~ (Aqcd/h) 2 times the ideal gas term, eq. (10). 
This is typical of a nonperturbative correction, as an inverse power of a (hard) 
mass scale. When /i ~ Aqcd, this is a large correction. When ji ^> Aqcd, 
numerically this is a very small contribution to the total pressure. Even so, as 
particles at the edge of the Fermi surface are the lightest excitations, even at 
large /i baryons dominate processes with low momenta. This implies that at 
large N c , phenomena involving the Fermi surface, such as superconductivity 
and superfluidity, are properly described by baryons, and not by quarks, for 
all \i ~ 1. 

By considering gluonic probes, it is clear that the theory is in a confined 
phase for T < Td and fi ~ 1. As discussed by Greensite and Halpern [4], 
at zero temperature the Wilson loop is insensitive to quarks at large N c : it 
exhibits an area law, with a nonzero string tension. Screening due to quarks 
enters through corrections ~ 1/N C . Adding a Fermi sea of quarks doesn't 
change this, as long as /i ~ 1. 

At nonzero temperature, the order parameter for deconfinement is the 
renormalized Polyakov loop [21,22]. Quarks induce an expectation value for 
the Polyakov loop, but this is ~ 1/AT C , versus a value ~ 1 in the deconfmed 
phase. Thus up to corrections ~ l/A^, as a function of temperature the 
expectation value of the renormalized Polyakov loop is independent of /i (for 
ji ~ 1); e.g., in fig. (1), T d ((i) is a straight line. 

That the theory confines can also be seen by exciting a quark, in the Fermi 
sea, with some external probe. If the quark is deep in the Fermi sea, knocking 
it out takes a probe with large momentum. When fi ^> Aqcd, a t first the 
resulting quark propagates like a hard quark. It, and the remaining hole in 
the Fermi sea, then scatter off of other quarks in the Fermi sea, knocking some 
out, creating other holes, and so on. Eventually, one ends up with not a single, 
unconfined quark, but a perturbed Fermi sea, characterized by some number 
of excited baryons and their holes. That the theory confines is clearer if the 
external probe carries momentum ~ Aqcd'- then one immediately sees that 
the only particles (and holes) excited near the Fermi surafce are not quarks, 
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but baryons. 

The above discussion applies to quarks of any mass, and leads to the phase 
diagram of fig. (1). An example is provided by the solution of QCD in 1 + 1 
dimensions [1]. In two dimensions there is only a confined phase, = oo. 
Schon and Thies showed that at nonzero quark density, the quark propagator 
remains infrared divergent, and so confined, as it is in vacuum [20]. This 
is exactly what one expects of a quarkyonic phase: that only color singlet 
excitations, such as mesons and baryons, have finite energy. It would be 
interesting to perform more detailed calculations, such as of the free energy, 
and how the properties of mesons and baryons change with [i. 

Returning to four dimensions, the crucial question is: where is the chiral 
phase transition for light quarks? For \x < m q , it surely coincides with the 
deconfining phase transition, and occurs at T d . We suggest, however, that the 
two transitions no longer coincide when \i > m q . If we take the Skyrme crystal 
as a guide to the quarkyonic phase, then chiral symmetry restoration occurs 
not at the mass threshold, but above \x = m q . This is because there must be 
some significant density of baryons in the Fermi sea to drive the transition. 
Further, the transition occurs when \i ~ m q , on the order of the constituent 
quark mass, and not at some fi which is asymptotically large in a power of 
N c . Schematically, this gives the blue line in fig. (1). 

This can also be seen in illustrative models. As an example of a possible 
solution to the Schwinger-Dyson equations [23], Wagenbrunn and Glozman 
[27] studied a model with a confining gluon propagator, ~ l/(k 2 ) 2 in mo- 
mentum space. At infinite N c , solutions to the Schwinger-Dyson equations 
are those where the quark propagator (and its vertex with gluons) change, 
but the gluon propagator doesn't. Computing with the quark propagator at 
fi 7^ 0, but leaving the gluon propagator unchanged, it is not difficult to see 
that increasing \x drives chiral symmetry restoration at some nonzero value of 
[i above the mass threshold. 

This can also be seen from the eigenvalue distribution of the Dirac operator. 
In vacuum, Banks and Casher showed that chiral symmetry breaking is driven 
by a nonzero density eigenvalues at zero eigenvalue [28]. When [i ^ 0, for 
Nf > 3 the eigenvalues spread out in the complex plane, and there is no 
simple analogous condition [15]. Even so, it is most natural then as \i increases 
above the mass threshold, that whatever effect the gauge fields have on the 
eigenvalues, that eventually it is overwhelmed by \l ^ 0. An explicit example 
of this is chiral symmetry restoration in random matrix models [15]. 

The splitting of the deconfining and chiral transitions can be represented 
naturally in effective models. Mocsy, Sannino, and Tuominen [7] showed that 
if the coupling between the Polyakov loop and the chiral condensate has one 
sign, the transitions coincide; for the other, they diverge. Thus this coupling 
vanishes at the point where the transitions diverge. 

At large N c , because of the changes in the magnitude of the pressure, the 
deconfining transition is expected to be of first order. Once the chiral transi- 



11 



tion no longer coincides with deconfinement, its order is presumably controlled 
by the usual renormalization group analysis [7], and depends strongly upon 
the number of flavors. 

There is one caveat which must be noted. Chiral symmetry, and its possible 
restoration, can be sensitive to pairing near the Fermi surface. Certainly for 
/i ^> Aqcd, one expects that quarks deep in the Fermi sea are best described 
as chirally symmetric. It is possible, however, that the baryons near the Fermi 
surface may experience non-perturbative effects which cause them to pair in a 
chirally asymmetric manner. This effect will manifestly be small, suppressed 
at least by ~ (A QC d/h) 2 - 

Up to this point, we have assumed that /i ~ 1, so that gluons are blind to 
quarks and their Fermi sea. Especially to understand finite N c , however, it is 
also necessary to consider values of \i which grow with (fractional) powers of 
N c . We assume that the temperature T ~ Aqcd, like 7^(0), the temperature 
for deconfinement at ji = 0. In perturbation theory, the pressure includes 
terms as 

PpertX^T) ~ N c N f /i 4 F , N c N f ^ 2 T 2 F 1 , N 2 C T A F 2 . (9) 

In the limit of large N c , Fq, F\ and F 2 are functions of the coupling constant 
g 2 and Nf. The coupling runs with both mass scales, fi and T. Thus when 
/i grows like a power of N c times Aqcd, perturbation theory in g 2 is a good 
approximation. As always, this is excepting power like corrections from the 
region near the Fermi surface. 

Consider /i ~ nI'^Aqcd- In this region, the quark contribution to the 
pressure, ~ N c fi A F ~ N 2 F , is as large as that of deconfined gluons, ~ 
N 2 F 2 . In this regime, the quark contribution to the pressure is independent 
of temperature, since ~ N C ^ 2 F\ ~ N^ 2 F\ is down by ~ l/\fW c . It is 
possible that there is a temperature dependent term in the pressure, induced 
by baryons near the Fermi surface; e.g., ~ N c NfT 2 Aq CD , which is down by 
~ 1/N C to the leading term, ~ N 2 . 

1 /2 

At larger values of fi ~ N c Aq C d, quarks contribute to the Debye mass in 
the limit of large N c , eq. (1). In this regime, the pressure is completely 
dominated by that of quarks at zero temperature, ~ N c {i 4 Fq ~ N% F . 
In perturbation theory, the gluon contribution to the pressure remains ~ 
N C T A F 2 ~ N 2 , with the temperature dependent part of the quark pressure 
also ~ N c fi 2 T 2 F 2 ~ N 2 ; thus both are down by 1/N C , relative to the quark 
term at zero temperature. 

1/2 

When /i ~ N c ' Aq C d, gluons are screened by quarks, and one is in a 
qualitatively new regime. At zero temperature, the Wilson loop no longer 
exhibits an area law; at nonzero temperature, the renormalized Polyakov loop 
acquires an expectation value of order one. Consequently, eventually the first 
order phase transition for deconfinement ends. It can do so in one of two 
ways: the phase boundary for deconfinement can either bend over to zero 
temperature, or it can end in a critical end point at T ^ 0. Our simple 
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arguments cannot predict which occurs. It is reasonable that either is only 

1 /2 

possible once gluons feel the quarks; i.e., when /x ~ N c ^qcd- If so, and 
there is a critical end point, the critical behavior is a correction in 1/N C to the 
total pressure, which is dominated by the zero temperature term for quarks. 

We remark that if Nf goes to infinity with N c , then all terms in eq. (9) 
are ~ N%. Indeed, when both iV c and Nf are large, then even in vacuum one 
cannot speak, rigorously, of confinement. We have no insight into this limit. 

Deryagin, Grigoriev, and Rubakov showed that at large N c , the color sin- 
glet pairing of chiral density waves dominates over the di-quark pairing of color 
superconductivity [12]. Their perturbative analysis is reliable for /x ^> Aqcd, 
as long as the quarks do not lie within ~ Aqcd of the Fermi surface. When 
they do, the pairing of baryons also contributes. 

Son and Shuster [12] showed that even for extremely large values of N c , 
Debye screening disfavors the pairing to chiral density waves, and quark color 
superconductivity dominates. In sec. (7), we adopt a similar criterion to esti- 
mate when in QCD there is a transition from a quarkyonic, to a perturbative, 
regime. 

6 Confinement and Chiral Symmetry Breaking 

At zero temperature and density, the effect of anomalies usually ensures that 
any confined phase is one in which chiral symmetry is broken. We now give a 
heuristic argument as to why this need not be true at nonzero density. 

Casher, and then Casher and Banks, argued that in the vacuum, confine- 
ment automatically implies the breaking of chiral symmetry [28]. Consider 
a meson, in which the quark propagates to the right, with a spin along its 
direction of motion. To remain a meson at rest, this must mix with a quark 
propagating to the left, which can happen by scattering off of a gluon. Its 
spin, however, is now opposite to the direction of motion, so its helicity has 
been flipped. Since in QCD the interactions preserve chirality, which for a 
massless field equals helicity, this change of direction cannot occur. It can if 
there is a mass condensate in the vacuum, which the quark can scatter off of, 
and flip its helicity. Note that this argument is especially tight in the limit of 
large N c , where the number of quarks in a meson is fixed. 

Now consider the similar process at nonzero temperature. Then besides 
scattering off of a gluon, one can scatter of a quark in a thermal distribution. 
However, if we consider the processes of both emission and adsorption, the 
total is n(Ek) — n^Ey); E% is the energy of the rightgoing quark, is Ey is 
the energy of the left going quark, and n(E) is the Fermi-Dirac statistical 
distribution function. 

For an isotropic distribution, as in thermal equilibrium, if the momenta 
are the same, then the two distribution functions cancel. Thus the process is 
only allowed when k ^ k! . Since the momenta of the right and left moving 
quarks are different, however, we end up with an excited meson, different from 
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the initial meson. That is, this process represents not a meson at rest, but 
scattering between a meson, and some thermally excited state, such as another 
meson. 

In general, this is fine in a thermal distribution: what we mean by a 
"meson" is a sum over states anyway. However, there is a problem at large N c : 
if T ^ and \x = 0, then all interactions vanish at large N c , and this process 
must be suppressed by powers of 1/N C . Thus Casher's argument suggests that 
at nonzero temperature, and [i = 0, the connection between chiral symmetry 
and confinement remains. 

This connection could be lost in the presence of a Fermi sea, however. The 
argument goes through as before, except now we scatter off of a quark in the 
Fermi sea. Physically, the quark in the test meson scatters off a quark in a 
baryon, which then scatters into a baryon hole. There is no inconsistency with 
the large N c expansion, because the scattering amplitude is large, of order one. 
This suggests that it is possible to have a confined, but chirally symmetric, 
phase at ji ^ 0. 

What of the constraints from anomalies, which after all, are due to ultra- 
violet effects? Certainly the anomaly itself is unchanged by temperature or 
density. However, their implications are less obvious when Tor/i are nonzero. 
Because of the breaking of Lorentz invariance at T and /i ^ 0, Itoyama and 
Mueller showed that many more amplitudes arise [28]. The anomaly relates 
these amplitudes, but not as directly as in vacuum. For example, Pisarski, 
Trueman, and Tytgat showed that the Sutherland- Veltman theorem, which 
relates the amplitudes for 7r° — > 77, does not apply at nonzero temperature 
[28]. Thus the connection between chiral symmetry breaking, and confine- 
ment, need not remain at nonzero T or fi. 

The above generalization of Casher's argument suggests that a system 
at fi 7^ is uniquely different from /1 = and T ^ 0. This may arise as 
follows. Anomalies are saturated by excitations with arbitrarily low energies; 
for example, when chiral symmetry breaking occurs, by pions. To model a 
chirally symmetric phase, consider massive, parity doubled baryons [27]. In a 
thermal distribution with [i = 0, massive modes are Boltzmann suppressed, 
and cannot be excited at low energy. At nonzero density, however, a Fermi sea 
of massive particles can be excited, with arbitrarily small energy, by forming 
a particle hole pair. 

We conclude this section by noting that the Skyrme model provides a direct 
example of a confined theory which satisfies the anomaly conditions. There, 
the Wess-Zumino-Witten term automatically incorporates all effects of the 
anomaly, such as 7r° — > 77, etc. This happens whether the background field is 
chirally asymmetric, as for large R, or chirally symmetric, as for small R. In 
either case, fluctuations from the Wess-Zumino-Witten terms automatically 
incorporate all anomalous amplitudes. 
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7 Quarkyonic Matter in QCD 



The large N c limit we consider is only simple because Nf is held fixed as 
N c — > oo. Since N c = Nf = 3 in QCD, then especially at fi ^ 0, it is far from 
clear that QCD really is close to this limit of large N c , and small Nf. For the 
purposes of discussion, we henceforth assume that it is. 

To extend the analysis at large N c to QCD, the principal physical effect 
to include is that of Debye screening. From (1), the Debye mass is rn 2 Debye = 
(2Nf/ir)a s (fi)fi 2 , where a s = g 2 /(4ir). This is to be compared with the scale 
of confinement, Aqqd- The latter is only approximate: probably a better 
measure of the confinement scale is not Aqqd V er se i but the mass of the p 
meson, ps 1 GeV. 

The Debye mass can be computed, at zero temperature and nonzero den- 
sity, to higher order in perturbation theory [6] . The really essential question is 
to know how the effective coupling runs. At nonzero temperature, and \x = 0, 
Braaten and Nieto suggested in the imaginary time formalism, as energies 
are always multiples of 2ttT, perhaps the effective coupling runs in the same 
way [24]. This was confirmed by computations to two loop order [24]. This 
implies that while ~ 200 MeV is relatively low, that the effective coupling 
is moderate in strength, even down to T d [21]. 

There doesn't appear to be any similar factor at nonzero density; at T = 0, 
the coupling should run like a s (fi/ (cAq C d)), where c is a number of order one. 
For purposes of discussion, we that assume perturbation theory is reliable for 
/i > 1 GeV; at this scale, the Debye mass is also ~ 1 GeV. 

A Fermi sea first forms when the quark chemical potential /i > M^/3 ~ 
313 MeV. Large N c suggests that dilute baryons persist only in a narrow win- 
dow, ~ Aqcd/N 2 . Then, at some scale above this, QCD becomes quarkyonic, 
in that the pressure rises rapidly. Notice that the increase in pressure is not 
associated with a phase transition. In terms of baryons, it appears to be due 
entirely to their strong interactions. Below fi ~ 1 GeV, it can also be viewed 
as due to the strong interactions amongst highly non-ideal quarks. 

Once Debye screening becomes significant above /i ~ 1 GeV, gluons are 
shielded, and the coupling may be (relatively) moderate in strength. Be- 
cause of Debye screening, at large \i scattering within Aq C d of the Fermi sur- 
face should be under control. This is unlike large N c , where Debye screening 
doesn't contribute until /i ~ \/N~ c . The transition from a quarkyonic regime, 
to one which is perturbative in quarks and gluons, is presumably smooth, as 
there is no order parameter to distinguish one from the other. (Assuming that 
the deconfining transition doesn't extend down to T = 0, see below.) 

Needless to say, our estimates for /x in the quarkyonic phase are extremely 
crude. Understanding the lower bound on fi requires matching onto models 
of nuclear matter; perhaps it might help by matching onto models consistent 
with large N c counting. The upper limit can be pinned down better through 
higher order calculations in perturbation theory at ji ^ and T = [6]. 
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Fig. 2. Possible phase diagram for QCD in the plane of temperature and baryon 
chemical potential. The blue line in the quarkyonic phase indicates the chiral phase 
transition. There is a critical end point for deconfinement. 

A possible phase diagram is drawn in fig. (2); following phenomenology 
[14], we plot this ClS CIS db function of the temperature and the baryon chemical 
potential, /xg. If large N c is a reasonable guide to N c = 3, this should look 
something like fig. (1), except that the sharp edges are smoothed out. For 
example, below the mass threshold, Ta should change little with /i; this appears 
to be true from numerical simulations on the lattice [18]. Similarly, at large 
N c nuclear matter rapidly goes from a dilute phase, to one which is dense and 
quarkyonic. We indicate this in the figure by drawing the quarkyonic phase 
slightly above Mjv, the nucleon mass. 

We expect that the chiral phase transition occurs in the quarkyonic phase, 
well above the mass threshold. For QCD, at present numerical simulations 
on the lattice indicate that for small /i, the deconfining and chiral transitions 
coincide, and are crossover. A chiral critical end point may exist in the plane 
of T and \ib [13]. One might conjecture that if such a critical end point 
exists, that the deconfining and chiral transitions split from one another at 
that point. 

Speculating in this manner, in the quarkyonic phase, the latent heat asso- 
ciated with the chiral transition might be relatively small. Certainly at large 
N c , the large increase in pressure, ~ N c , is not tied to the chiral transition. 
The behavior of the chiral transition is very sensitive to the number of flavors, 
and possible restoration of the axial U(l) symmetry, though. 

Consider the deconfining phase transition, after it splits from the chiral 
transition. At fixed /i, as T increases, one goes from a confined phase of 
parity doubled baryons, to one of quarks and gluons. Deconfinement could 
either remain crossover, or perhaps become first order again (from the splitting 
point?). If it does turn first order, it will then have to end in a critical end 
point, now for deconfinement. Alternately, a first order deconfining transition 
could perist down to zero temperature. We indicate this uncertainty by the 
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question mark in fig. (2). 

How can the quarkyonic phase be studied? For the total pressure, one 
should use a description not in terms of baryons, but in terms of quarks. Ad- 
mittedly, they are highly non-ideal quarks, but there are hints to their possible 
behavior from numerical simulations, on the lattice, at nonzero temperature. 
At T ^ and [i = 0, the pressure can be characterized by a generalized (or 
"fuzzy") bag model. This is a power series in 1/T 2 times the ideal gas term 
[16,17,21]. At T = 0, and nonzero fi, this suggests 

-Pquarkyonic (A 4 ) = /pert ~ f^l f 2 ~ B + . . . (10) 

Perturbative corrections are subsumed into f peT t- Nonperturbative corrections, 
such as due to confinement, are included in /i 2 and B. Because of the term 
~ /ic/U 2 , the constant B need not agree with the usual MIT bag term, even 
in sign. We note that such a parametrization arises naturally from a Skyrme 
crystal. In the simplest model, what is equivalent to a conformally symmetric 
term ~ /x 4 arises from the Skyrme term, ~ k. Power like corrections then 
arise from the usual sigma Lagrangian, /i 2 ~ f%, etc. The pressure from 
this generalized bag model should then match smoothly onto that of nuclear 
matter, with no phase transition between the two. 

In contrast, in order to compute properties near the Fermi surface, it is 
necessary to consider effective theories of baryons. In a phase with chiral 
symmetry breaking, at low density these must match onto models of nuclear 
matter. In a chirally symmetric phase, the baryons are parity doubled. One 
possibility is to use Nambu-Jona-Lasino (NJL) models, not of quarks [10,11], 
but of baryons. Linear models of parity doubled baryons may also be of use 
[27]. Phenomenon such as superfluidity and superconductivity, and trans- 
port properties in general, are dominated by these states. For parity doubled 
baryons, the patterns of baryonic superfluidity and superconductivity will be 
significantly constrained by anomaly conditions. One might guess that the 
scales of baryon pairing in the quarkyonic phase is on the order of those in 
ordinary nuclear matter; i.e., that the gaps are small, tens of MeV. 

Schafer and Wilczek [9] noted that for three light flavors, there is conti- 
nuity between a nucleonic phase and one with quark color superconductivity. 
While chiral symmetry breaking is large in a nucleonic phase, it is also gen- 
erated by color-flavor locking. This suggests that for quarkyonic matter with 
three flavors and three colors, that one possibility is for the massive, parity 
doubled baryons to form (small) gaps which spontaneously break chiral sym- 
metry. This is the simplest way by which massive baryons, which are now 
only approximately parity doubled, can satisfy anomaly constraints at jj, ^ 0, 
although we suspect there are others. 

One way of computing the properties of a quarkyonic phase is to use ap- 
proximate solutions of Schwinger- Dyson equations [23]. These are, almost 
uniquely, the one approximation scheme which includes both confinement and 
chiral symmetry breaking. They do have features reminiscent of large N c : at 
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low momentum, if chiral symmetry breaking occurs, the gluon propagator for 
Nf = 3 is numerically close to that for Nf = 0. At present, solutions at fi 7^ 
assume a Fermi surface dominated by quarks; if quark screening is not too 
large at moderate /i, these models should exhibit a quarkyonic phase. 

On the lattice, it is well known that while gauge theories have a sign prob- 
lem at nonzero quark density when N c > 3, that numerical simulations can be 
done for two colors. Recently, these were done at T = and ji ^ for heavy 
quarks: they exhibit superfluidity at the mass threshold and deconfinement 
well above it [19]. These simulations could be extended to light quarks, to see 
if the phase diagram is anything like that of fig. (2): e.g., at low temperature, 
is chiral symmetry restored before deconfinement? 

While our analysis is crude, existing prejudice has been that the phase 
transitions for deconfinement and chiral symmetry are inexorably linked to- 
gether. Large N c suggests that by moving out in chemical potential, that 
potentially one has the chance to see the two transitions separate. This could 
happen at rather high temperature, near that for the deconfining transition 
temperature at zero density, and relatively low density, less than that for nu- 
clear matter. Experimentally, it is possible to move out in the plane in /i, at 
high T, by going to "low" energies, such as at critRHIC and FAIR. Thus these 
facilities may explore not just a chiral critical end point [13], but quarkyonic 
phases, including one which is confined, yet chirally symmetric. 
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